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Abstract 

We investigate whether enhanced gravitational scattering on smaU scales (< 0.1mm), 
which becomes possible in models with large extra dimensions, can establish statis- 
tical equilibrium between different particle species in the early Universe. We calcu- 
late the classical relativistic energy transfer rate for two species with a large ratio 
between their masses for a general elastic scattering cross section. Although the 
classical calculation suggests that ultra-light WIMPs (e.g., axions) can be thermal- 
ized by gravitational scattering, such interactions are considerably less efficient once 
quantum effects are taken into account on scales below the Compton wavelength. 
However the energy transfer rate in models with several extra dimensions may still 
be sensitive to trans-Planckian physics. 



1 Introduction 

Higher dimensional models suggest a phenomenologically interesting solution to the hierarchy 
problem. A higher dimensional Planck mass M/j just above the electroweak scale is in fact com- 
patible with the observed weakness of four-dimensional gravity {Mp ~ lO^^GeV) if we consider 
relatively large compactification radii. In braneworld models where gravity freely propagates in 
A + d dimensions L and the d internal dimensions share the same compactification radius i?, 
the latter is related to as 



Within length scales of order R gravity is genuinely 4 + d dimensional and the gravitational 
potential starts growing as ~ l/r^^'^. This rules out Mp, ~ 1 TeV for d = 1 models, while 
leaving a small window for d = 2 2 . A fundamental Planck mass of ~ TeV opens the intriguing 
possibility of recording quantum gravity signals at the next generation of accelerators, where 
those energy scales will be actually probed 

In this note we investigate whether stronger gravity on sub-mm scales could affect the thermal 
history of the early Universe due to the enhanced gravitational scattering cross-section predicted 
by these models, even at energies and temperatures much less than the electroweak scale. Below 
the quantum gravity regime, we consider relativistic scattering processes between species of 
very different masses m ^ M in the limit of small scattering angle/low momentum transfer and 
estimate the efficiency of such processes for establishing statistical equilibrium between species. 
More specifically, we will consider the heavy species as relativistic and in thermodynamical 
equilibrium with an average energy per particle ~ T and the light species initially decoupled 
and "cold" i.e. with an energy per particle <C T. This is the case of the QCD axions, which 
acquire an effective mass m ~ W~^eV at temperatures as high as a GeV. 

Ultra-light weakly interacting massive particles (WIMPs) such as the axion are candidates 
for the cold dark matter (CDM) that is responsible for around 25% of the energy density of the 
Universe today and play a central role in cosmological structure formation. A crucial requirement 
for these ultra-light WIMPs is that they are non-thermal. In particular the QCD axions [3] arise 
from coherent oscillations about the minimum of its potential of the spatially homogeneous axion 
field that begin when the temperature drops below about 1 GeV at the QCD phase transition. 
They are supposed to be effectively non-interacting at this temperature and thus remain at 
rest until local overdensities undergo gravitational collapse, beginning the process of structure 
formation, from the bottom up. Today axions are supposed to be present in the virialised dark 
halo surrounding visible stars in galaxies, with velocity dispersion of order 100 km s~^. 

The gravitational cross-section of ultra- light WIMPs, enhanced by the presence of large 
extra dimensions, has recently been proposed |5j as a natural mechanism for the CDM "self- 
interaction" , advocated in [0] to resolve the "cusp problem" jlj in the centre of galaxies (other 
proposed solutions are found e.g. in In this paper we consider another effect of the 

enhanced gravitational scattering. We consider whether gravitational scattering between ultra- 
light WIMPs and much heavier but relativistic species in the early universe (such as electrons 
or neutrinos) could establish a statistical equilibrium between different species at energies well 
below the fundamental Planck scale. If the energy transfer is sufficient to make the axions highly 
relativistic, then they would no longer be viable candidates for the CDM in models with large 
extra dimensions. 
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2 Classical gravitational scattering in higher dimensions on a 
brane 



In the rest frame of the heavy particle, the hght particle feels a central attractive gravitational 
force 

|FWI = ^^^|5^, M-<.<i?, (2) 

where v is the speed of the light particle in this frame, and 7 = (1 — the usual 

relativistic Lorentz factor. Note that we use this equation as our definition of Md. Equation @ 
is derived from the 4 + d dimensional propagator and the lower bound r > is the limit of 
validity of such a tree-level calculation. Above the compactification radius, r > R, on the other 
hand, the usual Newtonian inverse-square law is restored (or r > iAdS in the Randall- Sundrum 
model [S] with an Anti-de Sitter bulk if i? > iAds)- 

In the small scattering angle limit the transverse momentum pj- acquired by the light particle 
can be estimated as 

26 , |F| 7 Mm 







where b is the impact parameter. It will be useful to define the velocity-independent dimension- 
less parameter 

vpr M 



(4) 



that, in the relativistic limit ~ 1, is just the transverse momentum in units of the energy of the 
particle. As long as we consider impact parameters bigger than the Planck scale cut-off M^^, 
we have e < M/Md <C 1. 

In Appendix A we calculate the energy transfered to ultra-light WIMPs due to gravitational 
scattering with much heavier, but relativistic particles such as electrons with Lorentz factors 
7 ^ 1 in the cosmological reference frame. We find, from Eq. (|22|). an energy transfer rate 



^ ^ 7rn,i?[7'-l] / dhbe\ (5) 

Jbuv 



where n.^ is the number density of the heavy particles and e(6) is given by Eq. The en- 
ergy transfer rate in higher-dimensional gravity [d > in Eq. @] is dominated by scattering 
events with small impact parameters due to the steep rise in the gravitational force on small 
scales r <C -R. Thus the energy transfer is not sensitive to the IR cut-off, biR (here given by 
the compactification scale ~ R), but is sensitive to the UV cut-off, buv- Note that in 4 di- 
mensional theory [d = 0) the energy transfer depends on the cut-off scales only logarithmically 
log(6/R/6(7y). 

We must cut-off our perturbative calculation at least at the fundamental Planck scale, buv ~ 
M^^, which is the limit down to which we can trust Eq.©. Note that for elementary particles 
with M < Md, the Schwarzchild radius rs ~ M^^(M/M£))^/^'^+"^^ of the heavy particle is much 
smaller than . A more serious problem is that our result is sensitive to trans-Planckian 
scattering and a reliable calculation requires knowledge of non-perturbative quantum gravity on 
smaller scales! The classical energy transfer rate due to gravitational scattering in brane-world 
models is sensitive to trans-Planckian physics. 
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On the other hand the Compton wavelength of our particles are much larger than so 
a full calculation requires a quantum mechanical treatment. We estimate the effect of quantum 
corrections in the following section, but in the rest of this section consider the effect of classical 
scattering. 

We can estimate the classical rate of energy transfer to the ultra-light particles as 

where we take the heavy species to be relativistic {v ~ 1) and in thermodynamical equilibrium 
(7M ~ T). Note that this result is independent of the rest mass of the heavy particles and 
depends only upon the number density of relativistic particles. 

Thus we would expect a rapid rise in the energy of ultra-light WIMPS due to classical 
gravitational scattering at sufficiently early times when E/E S> H. Note that Eq. ® has been 
derived in the limit 7m <C M, which was used in Appendix A to assume negligible recoil of the 
heavy particles during each collision. In practice the heavy particles are not infinitely heavy, and 
hence are not an infinite energy source, and the energy transfer rate must decrease as statistical 
equilibrium is approached. 

From © the relaxation rate of our light species can be expressed as a function of the 
temperature by taking the number density = LSAZ/T^/vr^ of fermionic relativistic species. 



Mf being the effective number of relativistic fermionic species. Gravitational interactions can 
efficiently establish equilibrium if, at some epoch in the early Universe, E/E ^ H, where the 
Hubble rate in the radiation dominated era is given by 

F=.10-^«(^)'GeV. (8) 



A more detailed calculation (see Appendix B) shows that the light particles acquire enough 
energy so that they are still relativistic at the beginning of matter domination Teq ~ 10 eV, 
if E/E > aH, where a ~ 30. In this case even if the light particles are initially "cold" when 
produced at Ta, they soon become relativistic due to the gravitational scattering by heavy 
relativistic particles in thermal bath. They remain relativistic until matter-radiation equality if 

Thus, for Md < 10 TeV, axions produced at Ta ~ GeV would acquire sufficient energy that, 
even by the time of matter-radiation equality, they are still relativistic and would therefore be 
ruled out as dark matter candidates. ^ 



^Note that our assumption of elastic scattering of the ultra-hght particles in the rest frame of the heavier 
particles required 7m < M which corresponds to i? < /T. For instance considering gravitational scattering 
with electrons at T ~ 1 GeV our calculation requires E/T < 10^^ . This is just sufficient to show that axion with 
mass m 10"'^ eV remain relativistic at Tcq ~ 10 eV. In practice we don't expect energy transfer to shut-off 
completely for /T < E <T, but it may become less efficient. 
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The reason why gravitational scattering remains efficient far below the fundamental Planck 
scale is the hierarchy between the fundamental Planck scale and the effective four-dimensional 
Planck scale, Mp ~ 10^^ GeV, which determines (in practice, suppresses) the Hubble expansion 
rate at these energies. For models where the fundamental Planck scale is as low as a TeV, 
gravitational scattering of ultra-light WIMPS remains efficient down to temperatures as low as 
lOMeV. This is analogous to the neutrinos remain coupled to baryons down to energies of order 
1 MeV even though the electroweak scale is of order 100 GeV. 



3 Quantum effects 

In calculating the energy transfer in the previous section we considered the particles as point-like 
and the scattering as a completely classical processes. By taking bjjv in ® as we have im- 
plicitly assumed that the two particles can come arbitrarily closed to each other, neglecting their 
finite Compton wavelengths. In the mass frame of the heavy particle the Compton wavelength 
of the light particle is {'ym)~^. By taking 7 as given by H14|l and by averaging over the angle 
9 we can estimate the Compton wavelength of the light particle as ~ M/{ET). By taking this 
rather than the Planck scale as the small scale cut-off buy in ® we obtain the energy transfer 

dE ^2d+1^2d+5 

here we assume R > (7771)"^ so that we still probe the higher-dimensional gravity. The relaxation 
rate E/E for axions is now much smaller that the Hubble rate H for Mpi ~TeV: 

E/E ^ MpE'^T''^+^ ^ in7-6<^^^V'^7^V'"''m'' nn 

H - M^dM^^+^ ~ \GeV) \Md) \m) ' ^ ' 

The last equality holds at T = Ta when the energy of the axions is E m. 

Rather than introducing an abrupt cut-off at the Compton wavelength, we can try to esti- 
mate the suppression due to the interference of the wavefunction for impact parameters smaller 
than the Compton wavelength, /c, of the particle. In the case of a uniform potential well of 
radius a the quantum mechanical cross section is suppressed, with respect to the classical one, 
by a factor {a/lc)^ 1^, i.e. (Tquantum = o"ciassicai(«/^c)^- The cross section has a direct inter- 
pretation in terms of "number of particles scattered". Therefore, in the integral in when 
considering impact parameters smaller than the Compton wavelength, b < Ic, we should include 
the corresponding suppression factor: 



dE o , /■ fbTE\^ , „ E'^T 

— ~ TTn^E[Y 
dt 



1] [ (EE.) be'db ^ ^'^l' I b-^'+'db. (12) 



Note that, as recently emphasized also in for d ^ 3 the steepness of the gravitational 
potential wins over the quantum suppression effect, and the process is still sensitive to trans- 
Plan ckian physics. For d < 3, on the other hand, most of the contribution comes from impact 
parameters larger than the Compton wavelength. Note also that for d ^ 3 the energy rate is 
given by (fTU)) with d = 3. 
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4 Conclusion 



In this letter, we studied the impact of an enhanced gravitational scattering cross-section in the 
early universe in the context of the brane-world models with a low fundamental Planck scale. 
In the presence of a general interaction of some given elastic scattering differential cross section, 
we calculated the relativistic energy transfer rate for two species with a large ratio between 
their masses. Due to the enhanced gravitational scattering cross section, ultra-light WIMPs 
could reach a statistical equilibrium with heavier relativistic particles in the thermal bath even 
well below the fundamental Planck scale. For example, the QCD axions produced at ~ 1 GeV 
could soon become relativistic and they remain relativistic until matter-radiation equality due 
to classical gravitational scattering. Axions are usually supposed to be non-interacting because 
their interactions are suppressed due to the existence of a much larger mass scale fa ~ 10^^ GeV. 
But in models with large extra dimensions their gravitational interactions can only be suppressed 
by the fundamental Planck scale, <C fa- 

However the classical energy transfer is dominated by scattering events with impact param- 
eters much less than the Compton wavelength of the light particles. If we cut-off our classical 
calculation at the Compton wavelength the effect disappears completely . If instead we attempt 
to model the quantum suppression of the calculation below the Compton wavelength we find 
that for d > 3 the result remains sensitive to the ultra-violet cut-off of the theory at the funda- 
mental Planck scale This is due to the steep rise in the gravitational force on small scales 
in higher dimensional spacetime. Thus, in a higher dimensional theory, the energy transfer due 
to gravitational scattering is sensitive to trans-Planckian physics. 

It is interesting to consider whether other phenomena could be sensitive to enhanced grav- 
itational scattering. We commonly neglect the effect of gravitational scattering as we expect 
it to be weak compared to all other interactions. But in models with large extra dimensions 
gravity can be much stronger on small scales than we naively imagine. And in contrast to four- 
dimensional gravity, the scattering is sensitive to trans-Planckian physics. Even the absence 
of detectable gravitational interactions could place constraints upon gravitational scattering on 
sub-Planck scales. 

Acknowledgments 

It is a pleasure to thank Paolo Nason, Gabriele Veneziano and Alberto Zaffaroni for useful 
discussions and correspondence. 

A Energy transfer rate by classical relativistic scattering 

We want to study the energy transfer between two particle species, with masses mi ^ m2, in 
terms of classical relativistic scattering. We assume that in the cosmological reference frame (e.g. 
the rest frame picked out by the cosmic microwave background) these species have speeds vi 
and V2, and that in this cosmic rest frame the velocity distribution of both species are isotropic. 

A given scattering event is characterized by the angle 9 between the initial 3-velocities vi 
and V2 of the two particles considered. Without loss of generality we can write the initial 
four-momenta in the cosmic rest frame as mi7i(l, vi, 0, 0) and 771272(1, V2cos0, V2sm6, 0) 
respectively, where 71 and 72 are the usual relativistic Lorentz factors. 
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In the particle-1 rest frame (RFl), the four-momentum of particle 2 is 
Pim = {E, Px, Py, Pz) = m272 (7i[l-^i^^2Cos6l], 71 [-^i + ^2 COS 61] , V2sme, 0), (13) 
from which the Lorentz factor 7 and relative velocity v in the RFl frame are derived as 

'^(^) = /, "^-.^so = 7172(1 - vif2Cos6'). (14) 

In RFl we can characterise the geometry of each scattering event by the two parameters b 
and <I>. We assume that we know the scattering deflection angle Q{b,v) in RFl as a function 
of the impact parameter b and velocity v. (In the case of interest in this paper ~ e/v^ ^ 1 
and e is given in equation @ as a function of b). After scattering, the light particle can acquire 
an orthogonal velocity component. We define an angle $ as the angle between the plane x — y 
and the plane containing 3- velocity of the light particle after scattering as well as particle 1. 
Only for <I> = or <I> = vr, the particle is scattered on the x — y plane and in general its velocity 
acquires an orthogonal component proportional to sin$. The final four-momentum after an 
elastic collision in RFl is thus given by 



Pfin = yE, px cos Q + Py sin cos <I>, py cos @ — px sin G cos ypi + p^ sin Q sin $j (15) 

Because we are working in the limit 1712/1^1 — > we require that the final speed of particle 2 
is the same as its initial speed in RFl and only its direction is changed. Total momentum is 
conserved due to the recoil of particle 1, but the kinetic energy transfered in RFl is negligible 
as m2/mi 0. 

For each particle of type 1 the differential rate of such an event is 

dfi{b,9,<^>) = bv{9)dn2{9)dbd^ , (16) 

where dn2 is the RFl-number density of those type 2 particles whose direction in the CMB 
reference frame is comprised between 9 and 9 + d9. 

By Lorentz-boosting back to the cosmic rest frame we obtain the total energy of the light 
particle after such an event: 



E(b,9,^) = m27i 



7 + ( — — 7 ) cos + f If 272 sin 9 cos $ sin ( 

71 



(17) 



Thus the energy acquired after the scattering, AE = E — 77172, in the limit of small scattering 
angle, G ~ e/u^, is 



AE{b,e,^) = m27i 



<bf (12 \ e{b) ^ . ^ 
7 + -//i\9 ^1^272 cos <P sm 9 



(18) 



2v{9Y V71 7 ^{oy 

The differential scattering rate in the cosmic rest frame is also easily worked out from 

dT2{9,b,^) = bv{9)dni{9)dbd^ . (19) 
By the assumption of isotropy in the cosmic rest frame 

dni{9) = ^-siii9d9, (20) 
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where ni is just the number density of species 1. 

Finally, the energy transfer rate in the cosmic frame 

dE 
lit 



AEdT (21) 



can be calculated. By integrating over <I>, from to 27r, the ^-dependent term in H18|) averages 
to zero. The integration over 6, from to vr can be worked out in the ultra-relativistic limit 
7i 3> 1. We are left with 



dE 

Ik 



f^IR 

7rni72m2[7?-l + 0(7r^)] / dbbe^ . (22) 

J brrv 



In the case of interest in this paper we must impose a UV cut-off buy ~ below which 

the classical scattering amplitude will not be valid, and at long wavelengths there is a cut-off 
bjR ~ R beyond which we recover Newtonian gravity. 



B Energy evolution of the light species 

We now want to follow the Energy evolution of the axions down from the temperature Ta ~ 
GeV where they are given a mass m ~ 10^^ eV and they are basically at rest in the reference 
frame of the cosmological observers. The axions' energy is controlled by two main processes: 
the gravitational interaction with some relativistic species (electrons, positrons or neutrinos) 
described in the text and the energy redshift due to the expansion of the Universe: 



fATfMp. E 



m 



dE = ^ . ET^ + (-dT) (23) 



The first term in the parenthesis comes directly from eq. ©, where temperature has been used 
instead of time (t ~ 10~^Mp/T^ during radiation domination) as independent variable and 
n^, = l.ST^/vr^ has been also used. The second and third terms represent the energy loss by 
redshift. In the range of interest between Ta and the beginning of matter domination at T ~ 10 
eV the third term is irrelevant. By posing E{Ta) = mn one finds the solution 

Ei;El = I_ e-C(T3-T3)^ (24) 
m Ta 

where 

IbirdM^ 2,-Kd \Md J 

The "heating up" process of the axions is very efficient at the beginning and E(T) reaches a 
maximum at a temperature which is a fraction of Ta- From then on, axions cool down by the 
effect of the redshift; in this regime E is basically linear in T: 

^^^e^^K T«Ta. (26) 
m Ta 

Requiring that the axion be non relativistic at Tgq = 10 eV implies requiring the above number 
to be of order one at that temperature. We obtain CT^ > 20, i.e. 
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